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ABSTRACT 

We have analyzed changes in the acoustic oscihation eigenfrequencies measured over the past 
7 years by the GONG, MDI and LOWL instruments. The observations span the period from 
1994 to 2001 that corresponds to half a solar cycle, from minimum to maximum solar activity. 

These data were inverted to look for a signature of the activity cycle on the solar stratification. 
A one-dimensional structure inversion was carried out to map the temporal variation of the radial 
distribution of the sound speed at the boundary between the radiative and convective zones. Such 
variation could indicate the presence of a toroidal magnetic field anchored in this region. 

We found no systematic variation with time of the stratification at the base of the convection 
zone. However we can set an upper limit to any fractional change of the sound speed at the level 
of 3 X lO"'^^ 



Subject headings: Sun: helioseismology, interior, magnetic fields 



1. Introduction 

Changes in the frequency of the solar p-mode 
oscillations have now been observed for more than 
a decade. Such changes affect both the central fre- 
quencies, Vni, and the frequency splittings, ^Vnim 
of low degree (Jimenez-Reyes et al. 2001), inter- 
mediate degree (Howe et al. 1999; Dziembowski 
et al. 2000) and very high degree modes (Ra- 
jaguru et al. 2001). A number of mechanisms 
have been proposed to explain these variations on 
frequency. Libbrecht & Woodard (1990) have ar- 
gued, on the basis of observations of intermedi- 
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ate degree modes, that the source of the pertur- 
bations must lie near the solar surface. Gough & 
Thompson (1988) and Paterno (1990) concluded 
that magnetic fields located near the base of the 
convection zone, with strengths significantly lower 
than 10^ G have no observable effect on p-mode 
frequencies. The stability analysis for magnetic 
fields by Morcno-Inscrtis et al. (1992) has shown 
that fields with strengths significantly larger than 
10^ G cannot be stored in this region. 

The p-mode frequency variations track rather 
well the changes of the activity strength of the 
solar cycle with time. It is thus plausible that 
these frequency shifts are due to variations of the 
mean magnetic field near the photosphere (Gough 
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& Thompson 1988; Goldreich et al. 1991). The 
influence of thin magnetic fibrils on the frequency 
shifts has been investigated by Bogdan & Zweibel 
(1985). Another possible cause for these frequency 
shifts is the presence of sunspots during solar ac- 
tivity (Braun et al. 1992). The dominant effect 
of sunspots on the propagation of acoustic waves 
is believed to be the dissipation of the acoustic 
energy and therefore, should decrease their am- 
plitudes and lifetimes (Hindman & Brown 1998; 
Rajaguru et al. 2001), while it has been observed 
that frequencies increase with increased magnetic 
activity. 

In the work presented here, we explore a mech- 
anism first suggested by Kuhn (1988), in which 
the sound speed perturbation associated to the 
observed changes of the photospheric latitudinal 
temperature distribution might be responsible for 
the frequency shifts seen during the solar cycle. 
Changes in the temperature distribution are them- 
selves due to the heat transport through the con- 
vection zone induced by the solar dynamo. 

Let us suppose that there is a magnetic field 
anchored below the base of the convection zone. 
To maintain pressure equilibrium, the gas pressure 
and thus the density inside the magnetized region 
must be lower than in its surroundings. The mag- 
netized fluid will thus experience a larger radia- 
tive heating and therefore the temperature at the 
top of the magnetized region will increase. This 
will also induce a change in the temperature gradi- 
ent that could be large enough to make the region 
above the magnetic field convectively unstable. In 
such a scenario, the base of the convection zone 
woiild locally drop, allowing the magnetized fluid 
to ascend by convective upflows, transporting ex- 
cess entropy to the photosphere (Kuhn & Stein 
1996). 

Numerical experiments by Kuhn & Stein (1996) 
have shown that entropy perturbations in the deep 
convection zone can produce strongly peaked tem- 
perature changes in regions below r = 1 that have 
a substantial acoustic signature (where r is the op- 
tical depth). They have also shown that the ther- 
mal perturbations that account for the solar acous- 
tic variability are consistent with the observed so- 
lar irradiance and luminosity changes that occur 
during the 11 year solar cycle. 

Luminosity changes, even if no larger than 
0.1%, must come from the release of energy 



stored somewhere in the solar interior and must 
be accompanied by a change in the solar radius. 
The ratio between relative luminosity and radius 
changes, hereafter VK, can help estimate the lo- 
cation of the region where this energy is stored 
(Gough 2000, and references therein). Theoret- 
ical calculations indicate that W increases when 
increasing the depth of the source of the varia- 
tions in luminosity. For instance W « 2 ■ 10~^ 
if the source is located in the outer layers of the 
convection zone, while W ~ 0.5 if the source is 
located in the solar core. Unfortunately there is a 
large scatter in the observed values of W . Indeed, 
recent measurements of W range from 0.021 as 
estimated by Dziembowski et al. (2001) and Antia 
et al. (2000a) to an upper limit of 0.08 derived by 
Emilio et al. (2000). 

Kuhn (2000) estimated that a 0.1% luminosity 
perturbation integrated over a solar cycle corre- 
sponds to about 10^^ erg. If this energy origi- 
nates in the tachocline and if the tachocline thick- 
ness is 0.05 i?0, the associated relative variation in 
sound speed at that depth would be on the order of 
5c/c w 10~^ or 10~^. Fractional changes in sound 
speed as small as 10"'' are easily accessible by he- 
lioseismic inversion techniques. Some attempts to 
find solar-cycle variations of the sound speed as- 
phericity and the latitude-averaged sound speed 
have been carried out using MDI and/or GONG 
data (Dziembowski et al. 2001; Howe et al. 1999; 
Antia et al. 2000b; Antia et al. 2001; Basu & Antia 
2001). However, none of them found any system- 
atic variation of the solar structure at the base of 
the convection zone that could be associated with 
the presence of a local toroidal magnetic field. 

We have extended the previous analysis to the 
latest data available, including LOWL data. Only 
common modes to all data sets were used, in an 
attempt to obtain comparable and significant re- 
sults to all instruments. This way, we can give a 
robust upper limit on the temporal variations of 
the solar internal stratification during the period 
1994-2001. 

2. Inversion Technique 

The inversion for solar structure, in particular 
sound speed c and density p, are commonly based 

on the linearization of the equations of stellar oscil- 
lations around a reference model (Gough & Koso- 
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vichev 1988; Dziembowski et al. 1990; Gough & 
Kosovichev 1990). The differences of the struc- 
tural profile between the actual sun and a model 
arc linearly related to differences between the ob- 
served frequencies and those calculated using that 
model. This relation is obtained using a varia- 
tional formulation for the frequencies of adiabatic 
oscillations. A general relation for frequency dif- 
ferences is given by 




where Svni art' the frequency differences between 
the actual sun and the model for the mode with ra- 
dial order n and degree /, and the correspond- 
ing relative error. The sensitivity functions, or 
kernels, if^'p(r) and Kpj^{r) are known functions, 
that relate the changes in frequency to the changes 
in the model. The hmctions dc/c and Sp/p are 
the unknown parameters to be inverted, i.e.: the 
relative difference in the sound speed and the den- 
sity respectively, and i?© is the solar radius. The 
term £~^F{vni) in Eq. 1 is introduced to take into 
account the so-called surface uncertainties; these 
include the dynamical effects of convection on the 
oscillation equations, as well as non- adiabatic pro- 
cesses in the near-surface layers (see Dziembowski 
et al. 1990, and references therein). 

Following standard procedures, we represent 
F{vni) as a Legendrc polynomial expansion. is 
the inertia of the mode, normalized by the inertia 
that a radial mode of the same frequency would 
have (for more details, see Gough & Thompson 
1991). 

If we take 5vni as the differences in the observed 
frequencies at two different epochs, rather than 
the differences in frequency between the actual sun 
and the model, 5c/ c and 5p/ p represent the vari- 
ation with time of the sun's internal structure, as 
long as our underlying theoretical model is very 
close to the actual sun. 

The inverse problem defined by Eq. 1 is well 
known to be an ill-posed problem (Thompson 
1995), whose solution is not unique. It can be 
solved using inversion methodologies that can be 
classified in two different techniques: the regu- 
larized least-squares methods (RLS, see Craig & 



Brown 1986) and the optimal localized average 
methods (OLA, Backus & Gilbert 1968). Both 
methods compute an estimate of the solution at 
a target location from a linear combinations of 
the observables, given a mesh of target locations. 
We have developed a variant of the RLS tech- 
nique, that we call the optimal mesh distribu- 
tion (OMD), that optimizes the mesh of target 
locations to avoid undesired high-frequency os- 
cillations of the solution. This optimization is 
achieved by computing a priori the spatial res- 
olution of the solution from the set of available 
observables and their uncertainties (Eff-Darwich 
& Perez- Hernandez 1997). The smoothing func- 
tion is itself defined also from the spatial reso- 
lution analysis and it is weighted differently for 
each radial point. This method ensures that the 
smoothing constraint is properly applied over the 
optimal mesh. 

3. Observational Data 

The observational data consist of mode frequen- 
cies computed from time series spanning different 
epochs and observed with different instruments. 
Namely, 57 sets based on 108-day-long time series 
derived from the GONG instruments (Leibacher 
et al. 1996) and spanning May 1995 to February 
2001; 27 sets based on 72-day-long time series de- 
rived from the MDl instrument (Schou 1999) and 
spanning May 1996 to November 2001; and 6 sets 
based on 1-year-long time series derived from the 
LOWL instrument (Jimenez- Reyes 2000; Tomczyk 
et al. 1995) and spanning 1994 to 1999. 

In order to use consistent data sets, only the 
modes common to all the sets for a given instru- 
ment were taken into account. As a consequence, 
the low degree modes {I < 13) present in some 
GONG data sets had to be rejected. Also this 
selection reduces the number of MDI and LOWL 
modes by 30% and 4% respectively. 

The MDI and LOWL sets were further reduced 
to only include the modes common to both instru- 
ments. This was not done with the GONG data 
set due to the small amount of common modes 
present. Finally, and again for consistency, we de- 
liberately restricted range of degrees we included 
to correspond to the highest degree available in 
the LOWL data set {i.e., I < 100). 

For each instrument and for each mode we com- 
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putcd the temporal frequency average. We sub- 
sequently subtracted the respective averaged fre- 
quencies from each set, leaving us with frequency 
changes with respect to this temporal average as 
a function of epoch. For the GONG and MDI 
sets, we also computed averages corresponding to 
1-year-long epochs. Such averaging reduced the 
scatter of the data while producing data sets com- 
parable to the LOWL sets. 

4. Results 

Figure 1 shows the relative change of the sound 
speed as a function of radius inferred from 1-year- 
long MDI, GONG and LOWL sets. These profiles 
show no significant changes at the level of a few 
times 10~^. The precision and resolution of the in- 
version is good enough to detect small variations 
of the stratification at the base of the convection 
zone. This is demonstrated in Fig. 2, where we 
show the sound speed profiles inferred from the 
1996 averaged MDI, GONG and LOWL data sets 
as well as sound speed profiles obtained by invert- 
ing the same mode sets, but after injecting fre- 
quency changes that result from a perturbation in 
the sound speed (as small as 3 x 10~^ and 3 x 10~^) 
between 0.68 and 0.70 i?©. This figure indicates 
that a perturbation of the sound speed at the base 
of the convection zone on the order of, or slightly 
smaller than 5 x 10^''', can be detected with the 
current precision resulting from 1-year-long time- 
series. Perturbations on the order of 10~^ fall in 
the noise level of our inversions. 

In an attempt to find temporal variations of the 
solar stratification at the base of the convection 
zone, we computed the mean value of ^c/c in the 
radial interval 0.69 < t/Rq < 0.72. This interval 
contains not only the base of the convection zone, 
r « 0.7133 i?0 (Basu & Antia 2001), but also the 
tachocline, r w 0.691 i?Q (Corbard et al. 1999), 
both closely related to the toroidal magnetic field 
responsible of the solar cycle. The resulting values 
are shown function of time in Fig. 3, for the 
inversions based on 1-year-long sets for all three 
instruments (GONG, MDI and LOWL), as well as 
on the GONG 108-day-long and MDI 72-day-long 
sets. 

Inversion profiles inferred from any linear inver- 
sion technique always correspond to the convolu- 
tion of the underlying solution by the resolution 



kernel (Thompson 1995). Therefore, even if the 
mode set used in a sequence of inversions remains 
identical, the resolution kernels will, at some level, 
change with time since the uncertainties change 
with time. Such variation could produce an ap- 
parent temporal behavior of the inferred profiles 
that docs not correspond to a real variation of the 
underlying true solution. To quantify this effect, 
we computed the averaging kernels at r = 0.69 Rq 
for all the inversions based on 1-ycar-long data 
sets. These were convolved with an artificial sound 
speed perturbation of the form: 

. r 3 X 10-5 for 0.67 < r/R^ < 0.71 
— = \ 5 X 10-5 for 0.91 < r/RQ < 0.93 
^ [ otherwise 

(2) 

where the sound speed perturbation centered at 
0.92 i?Q attempts to reproduce the results found 
by Dziembowski et al. (2000); Antia ct al. (2001), 
while a second perturbation at the base of the con- 
vection zone is also introduced. The convolutions 
g^o (t) are then calculated in the following way: 

q,^{t)= [ K{r,ro,t)^-^ dr (3) 

The resulting values of qr„{t), at = 0.69 7?©, 
relative to the average qav of all the convolu- 
tions calculated for every instrument, are shown 
in Fig. 4. This figure demonstrates that for the 
data from all three instruments the effect of the 
changes in the observed uncertainties on the in- 
verted profiles is negligible, corresponding to levels 
well below lO"''. 

The averaging kernel corresponding to the solu- 
tion at r = 0.69 Rq obtained from the 1997 MDI 
data set is shown in the right panel of Fig. 5. This 
is well located and indicates that our radial resolu- 
tion corresponds to 0.04 i?©. We should also point 
out that the averaging kernels have an important 
negative non-local contribution near the surface, 
a feature that affect at some level the solution at 
the base of the convection zone. The effect of this 
non-local component of the averaging kernels was 
quantified, in the case of MDI and GONG data, by 
performing inversions with data sets expanded to 
higher degrees, i.e., up to I = 150. The averaging 
kernel at r = 0.69 Rq obtained from the MDI 1997 
data shows a substantial reduction of the negative 
component located at r w 0.90 Rq when including 
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higher degree modes. But the temporal behavior 
of the stratification at the base of the convection 
zone do not significantly differ when the MDI data 
sets arc expanded from I < 100 to I < 150, as il- 
lustrated in the left panel of Fig. 5. 

Therefore, by including the data from all three 
instruments, and after assessing the effects of tem- 
poral changes of the resolution of the solutions, we 
can safely conclude that there is not significant 
systematic variations of the stratification at the 
base of the convection zone at the level of 3 x 10~^, 
and that this upper limit is constrained by the 
scatter present in the data. 

Vorontsov (2000) analyzed MDI data spanning 
from 1996 to 2000 and found systematic varia- 
tions of the radial solar stratification with time, 
expressed as relative changes of radius of 2 x 10~^. 
Our results obtained with MDI data are in good 
agreement with those found by Vorontsov (2000), 
in the sense that there is a systematic variation 
in the relative sound speed difference that is well 
correlated to the magnetic activity in the Sun. 
The maximum in solar activity corresponds to the 
maximum variation in sound speed. Since these 
results are not seen when analyzing LOWL and 
GONG data they should still be taken with some 
degree of scepticism. 
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Fig. 1. — Relative sound speed differences inferred using 1-year- long averaged GONG (panel a), LOWL 
(panel b) and MDI (panel c) data, shown for clarity only for the years 1997, 1998 and 1999. The solutions 
are shown only where the inversions are reliable. For clarity, the error bars are only shown for the inversion 
of the 1997 data; the error bars of the other results are almost identical. 
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Fig. 2. — Relative sound speed differences inferred from the 1996 averaged GONG (panel a), LOWL (panel 
b) and MDI (panel c) data sets. For each panel, the solid line shows the inversion similar to those presented 
in Fig. 1, whereas the other two profiles correspond to inverted profiles from the same mode set but after 
the injection of a frequency change corresponding to a small perturbation in the sound speed between 0.68 
and 0.70 i?0, of amplitude e = 3 x 10~^ and 3 x 10~^ respectively. For clarity, error bars have not been 
included, although they are very similar to those shown in Figure 1. 
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Fig. 3. — Temporal variations of the average of 6c/c at ro — 0.69 i?© computed over the radial points that 
lie within the interval 0.69 < t/Rq < 0.72. Left panel shows the results for LOWL, MDI and GONG 
1-year-long data sets. Right panels compare results for the 108-day-long and 1-year-long GONG sets and 
the 72-day-long and 1-year-long MDI sets. 
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Fig. 4. - Contribution of the time- varying averaging kernel to spurious sound speed profile changes with 
time; see text for details. 
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Fig. 5. — Left panel: temporal variations of the average of 5c/ c over the radial points that lie within the 
interval 0.69 < t/Rq < 0.72. Results correspond to 1-year-long MDI data with degrees expanding up 
to I = 100 (shown as diamonds) and I = 150 (shown as filled circles). Right panel: averaging kernels 
corresponding to the inversion of the 1996 data at r = 0.69 i?0. The solid line corresponds to the averaging 
kernel obtained from the inversion of data with degrees up to Z = 100, while the dotted lines represents the 
averaging kernel obtained from the inversion of data with degrees up to Z = 150. 
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